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摘 Xe. 针对 宁夏 中 部 干旱 带 压 砂 地 土壤 肥力 下 降 、 灰 绿 板 兰 矿 质 元 素 淋 溶 释放 规律 不 明 等 问题 ， 
通过 室内 模拟 灰 绿 板 兰 冻 融 干 温 循环 交替 试验 ,采用 修正 的 Elovich 方程 .抛物 线 方程 , 双 常 数 速 率 
方程 一 级 动力 学 方程 对 各 矿质 元 素 的 累积 释放 曲线 进行 拟 合 ,研究 不 同 粒 径 板 内 在 冻 融 干 温 特 
环 作用 下 的 释放 动力 学 特征 及 最 优 动力 学 方程 。 结 果 表 明 : 两 种 粒 径 灰 绿 板 崇 淋 溶液 中 矿质 元 素 


E 


淋 溶 总 量 随 循环 次 数 增 大 而 增 大 ,1 em ERR Kae AIK ARMIES BIR EA l dB 
环 次 数 下 均 大 于 3 cm 粒 径 。 灰 绿 板 岩 矿质 元 素 淋 溶 释 放 是 多 因素 共同 控制 的 物理 、 化 学 过 程 ,其 


释放 过 程 可 大 致 分 为 快速 反应 阶段 和 反应 趋 于 


FE 衡 阶段 ; 压 砂 地 土壤 中 Ca、K、Mg、P 元 素 的 供给 量 


运用 修正 的 Elovich 方 程 预测 结果 较 好 ,抛物 线 方程 更 适用 于 描述 S 元 素 的 释放 规律 。 研 究 结 果 可 
为 字 夏 中 部 干旱 带 旱 作 农田 土壤 肥力 调控 提供 决策 参考 。 
X 键 词 : 冻 融 干 湿 作 用 ; 灰 绿 板 岩 ; 矿质 元 素 ; 淋 溶 释放 
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灰 绿 板 岩 为 宁夏 中 部 干旱 带 压 砂 常 采用 的 砾 
石 ,将 其 覆 于 当地 砂 壤 士 上 可 对 土壤 起 到 蕾 水 保 
JH . 增 温 . 压 碱 等 作用 2 。 随 着 当地 种 植 结 构 的 改 
变 及 压 砂 年 限 的 增长 ,土壤 肥力 下 降 5、 砂 地 “ 老 
化 ”、 丢 弃 砂 田 面 积 增 大 “等 问题 日 益 凸 显 , 目 前 关 
于 压 砂 地 的 研究 大 多 集中 于 土壤 水 盐 调 控 \ 土 壤 结 
构 组 成 .有机质 含量 等 方面 ”” ,而 对 砂 地 表面 灰 绿 
板 岩 在 经 过 降水 、 灌 溉 等 多 年 的 冻 融 干 湿 往 复 循环 
作用 下 的 内 部 矿质 元 素 淋 溶 释放 特征 及 释放 动力 
学 模拟 研究 却 鲜 有 报道 。 

岩石 在 风化 过 程 中 其 元 素 的 淋 溶 释放 受理 化 
性 质 环境 等 因素 的 影响 ,岩石 处 在 不 同 环 境 中 元 
素 释放 均 有 所 差异 ""'""。 有 研究 表明 长 白山 森林 景 
WARE P Sr P, Ca, Na 元 素 的 淋 洲 过 程 较 明显 22 ,而 
在 南方 喀斯特 地 貌 地 区 Mg 等 活动 性 较 强 的 元 素 易 
发 生 淋浴 所 ; 灰 绿 板 岩 在 不 同 补 水 量 水 平 下 Mn Fe, 
Ba、Cl、Na、Si 元 素 淋 溶 量 存在 差异 ;煤层 围 岩 中 
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Cd 元 素 在 冻 融 条 件 下 的 释放 较 降雨 淋 溶 条 件 下 更 
容易 ,元 素 的 淋 溶 释放 受 粒 径 大 小 的 影响 "”。 还 有 
部 分 学 者 对 元 素 的 释放 动力 学 过 程 进行 了 研究 , 王 
继 纲 等 "9 采用 Elovich 方 程 和 Freundlich 修 正式 较 好 
的 描述 酸性 条 件 下 Cu Zn 的 动力 学 释放 过 程 ; 黎 晓 
起 等 "用 双 常 数 速 率 方 程 较 好 地 描述 有 机 酸 草 酸 
降解 海洋 玻 滩 物 中 重金 属 Cr 的 动力 学 过 程 。 综 上 ， 
关于 元 素 淋 溶 的 研究 大 多 集中 在 湿润 地 区 或 岩石 
长 期 处 于 湿润 状态 下 的 淋 溶 规律 , 极 少 关注 于 西北 
冻 融 干 湿 交 蔡 环境 下 的 元 素 淋 溶 规律 ;关于 元 素 淋 
溶 释 放 动 力学 的 研究 多 以 重金 属 元 泰 为 主 , 较 少 研 
究 矿 质 元 素 或 全 元 素 的 释放 动力 学 过 程 。 

岩石 长 期 处 在 富 含 0,.C0; 和 水 等 的 环境 中 易 
使 岩石 内 部 孔 际 增 大 、 矿 物 成 分 溶解 ,矿质 元 素 淋 
VER?" ,进而 直接 影响 土壤 的 营养 结构 。 灰 绿 
板 岩 在 风化 过 程 中 元 素 的 淋 溶 规律 及 其 淋 溶 释放 
过 程 是 否 影响 当地 土壤 肥力 .改变 土壤 结构 等 问题 
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尚 不 明确 。 鉴 于 此 ,本 研究 选择 压 砂 地 灰 绿 板 岩 为 
研究 对 象 ,通过 室内 模拟 冻 融 干 湿 循 环 试 验 ,探讨 
不 同 粒 径 .循环 次 数 对 灰 绿 板 岩 矿质 元 素 淋 溶 释放 
的 影响 并 建立 矿质 元 素 在 水 岩 作用 下 的 释放 动力 
学 模型 ,对 摸 清 压 砂 地 灰 绿 板 岩 元 素 的 迁移 淋 溶 规 
律 . 压 砂 地 生态 修复 和 可 持续 发 展 具 有 重要 意义 。 


1 材料 与 方法 


1.1 试验 设计 

试验 样品 选择 当地 压 砂 的 灰 绿 板 岩 。 聘 请 当 
地 村 民 在 雨水 冲击 河道 中 利用 钢 卷 尺 手 动 挑选 直 
径 为 1 cm、3 cm 的 灰 绿 板 岩 带 回 实验 室 ,并 用 蒸 馅 
水 冲洗 样品 去 除 表面 杂质 。 

本 研究 依据 《SL_T264-2020 水 利水 电工 程 岩 石 
试验 规程 CB/T 50226-2013 工程 岩 体 试验 方法 标 
准 》 设 置 粒 径 (1 cm、3 ecm) 和 循环 次 数 (20 次 .40 次 、 
60 次 80 次 ) 两 因素 室内 模拟 试验 ,试验 共计 8 个 处 
理 ,每 个 处 理 重复 3 次 ,共计 24 次 试验 。 冻 融 干 湿 
交替 条 件 采 用 烘箱 和 冰箱 控制 。 

根据 4SL_T264-2020 水 利水 电工 程 岩 石 试验 规 
程 ) 及 人 研究 区 夏秋 季 干 湿 和 冬 春季 冻 融 特征 , 干 湿 
状态 模拟 参考 研究 区 1989 一 2018 年 干 湿 期 (6 一 11 
月 ) 夏 秋季 累 年 月 平均 最 高 气温 ,并 结合 当地 气候 
条 件 设 置 干燥 温度 为 30 9C. ,湿润 温度 为 4 %C; 冻 融 
状态 模拟 依据 研究 区 1989 一 2018 年 冻 融 期 (12 月 一 
次 年 5 月 ) 冬 春季 累 年 月 平均 最 低 气 温 ,设置 冻结 温 
度 为 -20 % ,结合 1989 一 2018 年 冬 春季 累 年 月 平均 
最 高 气温 设置 溶解 温度 为 25 %C; 干 湿 冻 融 时 长 根据 
日 降雨 量 及 最 低 气 温 等 因素 综合 考虑 设置 为 干燥 
时 长 为 22 bh, 泡 水 溶解 时 长 为 2h, 冻 结 时 长 174, 溶 
解 时 长 为 7h; 一 个 干 湿 冻 融 循环 时 长 为 48 h ,进行 
80 次 循环 ,每 循环 20 次 测定 相应 指标 。 

1.2 指标 测定 

1.2.1 矿质 元 素 含 量 测定 ” 灰 绿 板 岩 淋 溶液 中 存在 
多 种 矿质 元 素 ,其 中 含量 最 多 的 元 素 为 Ca K, Mg, 
S.P, 这 些 元 素 是 植物 生长 所 需 大 量 元 素 ,对 植物 生 
长 具有 重要 影响 ”。 本 文 主要 对 这 5 种 元 素 的 淋 溶 
规律 释放 动力 学 过 程 进行 分 析 模 拟 。 

冻 融 干 湿 循 环 过 程 中 ,每 循环 20 次 后 取 50 mL 
灰 绿 板 岩 淋 溶液 样品 ,用 45 um 微 孔 滤 膜 过 滤 后 , 采 


JH ICP-OES 720 测 定 淋 溶液 中 矿质 元 素 的 含量 。 其 
中 ,矿质 元 素 淋 溶 总 量 为 每 循环 20 次 后 1L 淋 溶液 
中 Ca、K、Mg、S.P 元 素 的 含量 之 和 (mg*L')。 人 矿质 
元 素 累 积 淋 溶 量 为 矿质 元 素 淋 溶 总 量 随 循 环 次 数 
增加 的 累积 量 (mg:1L')。 
1.2.2 微观 结构 观测 “于 试验 前 后 采用 扫描 电镜 
(德国 蔡司 Zeiss EV0180) 观 察 灰 绿 板 岩 表面 微观 形 
SUERTE o 
13 数据 分 析 

(1) 各 矿质 元 素 在 快速 反应 阶段 的 释放 量变 化 
量 ( Ay ) ,计算 公式 如 下 : 

Ay = yay 7 Yo (1) 

式 中 : Ay Ay post x WE T E EK) EBLE ft mg L’); 
Yao 为 循环 40 次 时 的 元 素 累 积 淋 深 释 放量 (mg:I); 
yw 为 循环 20 次 时 的 元 素 累 积 淋 溶 释放 量 (mg*L'!)。 

(2) 本 文采 用 决定 系数 ( 尼 ) .平均 精度 (4 ) 对 方 
程 模拟 效果 进行 评估 ,计算 公式 如 下 : 


NS 


121 24070 (2) 
Y6.-» 
5-9 
A741 - M —.— > x 100% (3) 


StH: y ASME mg L); y, ATU (mg L’); 
y APP ISA (mg: L); NOJFEZS ZEE, ARK , ATE 
预测 的 精度 越 高 , 当 4>80% 时 ,说 明 方 程 的 预测 结 
果 可 反映 实际 情况 = 。 

(3) 参考 国内 外 学 者 的 研究 成 果 , 本 研究 拟 采 
用 动力 学 方程 ( 表 1)。 

(4) 使 用 Microsoft Excel 2010 软件 进行 数据 计 
算 与 分 析 ,Origin 2021b 拟 合 动力 学 方程 并 绘图 ， 


表 1 动力 学 方程 表达 式 


Tab.1 Expressions of kinetic equations 


(000 AER R HUBS 
一 级 动力 学 方程 Inyza- bx [17] 
修正 的 Elovich 方 程 y=atblnx [24] 

双 和 常数 速率 方程 Inyz a + blnx [25] 

(Freundlich 修正 式 ) 

抛物 线 方程 y=Q+Dx [26] 
iE: y 为 各 元 素 累 积 释放 量 (mg*1'); x 为 循环 次 数 ; a、5 为 常数 。 


下 同 。 
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PCAS 软 件 对 扫 面 电子 显微镜 图 片 进 行 分 析 处 理 。 


2 结果 与 分 析 


2.1 矿质 元 素 淋 溶 量变 化 规律 

由 各 循环 次 数 下 矿质 元 素 淋 溶 总 量变 化 规律 
(图 1) 可 知 ,两 种 粒 径 灰 绿 板 岩 麻 洲 液 中 矿质 元 素 
淋 溶 总 量 随 循环 次 数 增 大 而 增 大 。40 次 .60 次 、80 
次 循环 后 ,两 种 粒 径 灰 绿 板 岩 淋 溶液 中 矿质 元 素 淋 
HE EAE TE 0.001 水 平 存在 极 显著 差异 ,说 明 灰 绿 板 
宕 粒 径 的 大 小 在 一 定 程度 上 影响 其 矿质 元 素 淋 溶 
总 量 。 两 种 粒 径 灰 绿 板 岩 在 60~80 次 时 处 在 一 个 绥 
增 阶 段 但 尚未 达到 峰值 ,说明 冻 融 干 湿 循 环 至 80 次 
时 岩石 内 部 仍 有 少量 矿质 元 素 未 溶解 释放 。 


9000F 3 cm 粒 径 灰 绿 板 岩 
8000 上 1 cm 粒 径 灰 绿 板 岩 eK ane 


7000 上 


*** 


6000 
5000 


矿质 元 素 淋 溶 总 量 /mg.L 
> 
E 
S 


3000 
2000 
1000 
0 
20 40 60 80 
循环 次 数 /次 


注 :** 表 示 在 P<0.01 水 平 上 差异 极 显 著 ;*** 表 示 在 P<0.001 水平 上 
差异 极 显 著 。 
图 1 两 种 粒 径 灰 绿 板 岩 矿质 元 素 淋 溶 总 量变 化 规律 


Fig.1 Variation rule of total leaching amount of mineral 


elements in grey-green slate with two particle sizes 


由 两 种 粒 径 灰 绿 板 岩 矿 质 元 素 累 积 淋 溶 总 量 
变化 规律 (图 2) 可 知 ,1 cm 粒 径 灰 绿 板 岩 矿 质 元 素 
累积 淋 溶 总 量 在 不 同 循环 次 数 下 均 大 于 3 em BE, 
说 明 粒 径 越 小 越 容易 使 岩石 内 部 矿物 溶解 ,导致 元 
素 水 迁移 能 力 和 淋 失 率 增 大 。 累 积 淋 溶 量变 化 曲 
线 的 斜率 可 在 一 定 程 度 反映 元 素 的 淋 溶 速率 ,在 冻 
融 干 湿 循环 过 程 中 1 cm 粒 径 灰 绿 板 岩 的 淋 溶 速率 均 
大 于 3 cm 六 径 灰 绿 板 岩 ,1 cm 粒 径 灰 绿 板 岩 体积 较 
小 ,相同 循环 次 数 下 ,水 溶液 进入 1 cm BLAIR ER AR 
岩 内 部 孔 际 和 和 裂 际 的 体积 较 多 ,反之 ,3 cm 粒 径 较 
Jy 在 冻 融 循环 过 程 中 1 cm 粒 径 受 冻 胀 力 的 破 
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25000「 C 3 cm 粒 径 灰 绿 板 岩 

1 cm 粒 径 灰 绿 板 岩 
| 3 cm 粒 径 灰 绿 板 岩 变 化 趋势 pz 
~… 1 cm 粒 径 灰 绿 板 岩 变化 趋势 07 


15000 


10000 


矿质 元 素 累 积 淋 洲 总 量 /mg.L 


5000 


20 40 60 80 
循环 次 数 /次 
图 2 两 种 粒 径 灰 绿 板 岩 矿 质 元 素 累积 淋浴 总 量变 化 规律 


Fig.2 Variation rule of cumulative leaching amount of 


mineral elements in grey-green slate with two particle sizes 


坏 程 度 较 强 , 干 湿 循环 过 程 中 1 cm 岩石 内 部 在 水 的 
渗流 迁移 中 易 将 某 些 元 素 淋 出 。 
2.2 不 同 粒 径 板 岩 矿质 元 素 的 释放 动力 学 特征 

灰 绿 板 岩 在 冻 融 循环 干 湿 循 环 20 次 ,40 次 .60 
次 和 80 次 时 分 别 测定 矿质 元 素 释 放量 ,结果 见 图 
3。Ca KMg、.S 元 素 的 动力 学 曲线 形状 大 致 相 
同 ,循环 次 数 越 多 ,各 元 素 的 释放 量 越 大 。 两 种 粒 
径 灰 绿 板 岩 矿质 元 素 的 释放 过 程 可 分 为 2 个 阶段 : 
第 1 阶段 为 快速 反应 阶段 (循环 20~40 次 ) ,这 一 阶 
段 图 中 曲线 斜率 较 大 ,元素 释放 速率 较 快 ,在 干 湿 
冻 融 循环 初期 , 灰 绿 板 宕 表面 可 溶 和 易 深 元素 释 
放 , 使 淋 溶 液 中 元 素 含量 增 大 ;第 2 阶段 为 趋 于 稳定 
阶段 (循环 40~80 次 ) ,这 一 阶段 图 中 曲线 斜率 较 小 ， 
元 素 释 放 速 率 较 第 一 阶段 稍 缓慢 ,可 能 是 大 部 分 表 
面 易 淋 溶 元素 释 放逐 渐 减 小 ,元 素 从 岩石 内 部 的 微 
了 筷 际 中 缓慢 释放 ,导致 矿质 元 素 累 积 淋 溶 量 趋 于 

对 照 图 3 可知 1 cm 粒 径 灰 绿 板 岩 中 Ca、K、Mg、 
P、S 元 素 的 累积 淋 溶 释放 量 在 不 同 循环 次 数 下 均 高 
于 3 cm 粒 径 灰 绿 板 岩 的 累积 释放 量 。 在 快速 反应 
阶段 ,通过 计算 各 矿质 元 素 在 快速 反应 阶段 的 释放 
量变 化 量 ( Ay ) 可 知 ,1 em fife IK ARAL Ca. K Mg, 
PS 元 素 Ay 分 别 为 1172.315 mg: L!,2263.415 mg， 
L. 141.205 mg-L'',48.839 mg .二 .25.549 mg: L"'; 
3 cm ALTE JK EUR Ca, K Mg P,S7628 Ay 分 别 为 
370.226 mg-*L/,615.890 mg* L',360.054 mg: LL、 
18.427 mg*L/,13.697 mg-L, 由 以 上 数据 及 图 中 
曲线 的 斜率 可 得 出 1 em 粒 径 灰 绿 板 岩 的 各 矿质 元 
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图 3 灰 绿 板 岩 矿 质 元 素 释放 动力 学 曲线 


Fig. 3 Release kinetic curve of mineral elements in grey-green slate 


素 累 积 释放 曲线 的 斜率 相对 于 3 em 粒 径 灰 绿 板 岩 
的 斜率 更 大 ,说 明 初 始 阶 段 粒 径 越 小 灰 绿 板 岩 释放 
各 矿质 元 素 能 力 越 大 。 
23 不 同 粒 径 灰 绿 板 岩 矿质 元 素 的 释放 动力 学 方 
程 模拟 

内 外 学 者 对 元 素 释 放 动 力学 进行 研究 并 建 
立 与 之 相关 的 释放 模型 ”3 " 。 水 岩 作用 下 矿质 元 素 
的 淋 溶 释 放 是 一 个 复杂 的 物理 .化 学 等 相互 作用 的 
过 程 ,为 定量 描述 矿质 元 素 的 释放 动力 学 行为 ,分 
别 采用 一 级 动力 学 方程 .修正 的 Elovich 方 程 . 双 和 常 
数 速 率 方程 (Freundlich 修正 式 ) 抛物线 方程 对 矿 
质 元 素 累 积 释 放量 进行 拟 合 。 
2.3.1 一 级 动力 学 方程 ”各 矿质 元 素 的 释放 动力 学 
过 程 可 以 用 释放 量 的 对 数 与 循环 次 数 的 线性 关系 
来 表示 ,一 级 动力 学 方程 常用 于 描述 体系 中 扩散 机 
理 较 单一 的 过 程 , 比 如 颗粒 表面 的 扩散 等 "…。 采 用 
一 级 动力 学 方程 将 不 同 粒 径 矿 质 元 素 累 积 释放 量 


进行 拟 合 ,结果 见 表 2。 

由 表 2 可 知 ,3 em 粒 径 灰 绿 板 岩 各 矿质 元 素 累 
积 释 放量 采用 一 级 动力 学 方程 拟 合 的 尼 在 0.756-~ 
0.919 之 间 , 均 达到 显著 水 平 ;1 em BERE IK E Bo e 
矿质 元 素 累 积 释 放量 拟 合 的 尺 在 0.629~0.858 之 间 ， 
K PS 元 素 采 用 一 级 动力 学 方程 的 拟 合 达到 显著 水 
平 ,Ca Mg 元 素 不 显著 。 从 4 来 看 ,两 种 粒 径 灰 绿 板 
A FAA 1 cm 板 岩 中 S$ 元 素 的 4 高 于 80% ,其 余 均 未 
达到 80% ,说 明 该 方程 预测 不 同 矿质 元 素 的 实际 释 
2.3.2 修正 的 Elovich 方程 ”修正 的 Elovich 方 程 描 
述 体系 中 一 系列 反应 机 制 共 同 作 用 的 反应 过 程 ?， 
利用 修正 的 Elovich 方程 对 不 同 粒 径 矿质 元 素 的 释 
放 动 力学 过 程 拟 合 ,结果 见 表 3。 

从 屁 来 看 ,两 种 粒 径 灰 绿 板 岩 各 矿质 元 素 累 积 
释放 量 能 够 较 好 地 采用 修正 的 Elovich 方程 来 摘 
述 。3 cm 粒 径 灰 绿 板 岩 矿 质 元 素 累 积 释 放量 的 尼 


表 2 一 级 动力 学 方程 拟 合 


Tab.2 Fitting of first-order kinetic equation 


粒 径 /em 元 素 方程 决定 系数 (RR) 平均 精度 (4)/% 
3 Ca lny=5.259+0.017x 0.767 52.991 
K Iny=6.060+0.016x 0.811 64.396 
Mg lny=5.326+0.016x 0.756: 39.493 
P 1ny-2.52840.015x 0.768" 70.025 
S lny=2.418+0.018x 0.919" 75.061 
1 Ca Iny=6.459+0.015x 0.671 54.029 
K lny=7.274+0.014x 0.694 72.153 
Mg Iny=6.291+0.015x 0.629 62.362 
P lny=3.461+0.013x 0.661' 68.052 
S Iny=2.794+0.019x 0.858" 85.724 
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表 3 修正 的 Elovich 方程 拟 合 
Tab.3 Fitting of the modified Elovich equation 
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粒 径 /cm 元 素 方程 TRIER BUR’) 平均 精度 (4)/% 
3 Ca y=-1194.561+444.033Inx 0.978" 82.915 
K y7-1976.9724783.465l1nx 0.988" 86.779 
Mg y7-1029.610--398.6011nx 0.969 80.065 
P y=-52.029+20.999 nx 0.971 88.702 
S y7-62.433424.301Inx 0.973* 91.707 
1 Ca y7-2899.0321134.854lInx 0.923 90.219 
K y=-5543.136+2259.953Inx 0.932 94.999 
Mg y7-2681.82241025.8811nx 0.895 89.520 
P y=-107.518+45.157Inx 0.906" 90.292 
S y7-104.968--39.4101lnx 0.909* 91.695 


均 达 到 0.95 以 上 ,两 种 粒 径 矿 质 元 素 累积 释放 量 的 
尼 大 小 为 :K>Ca>S>P>Mg。 修 正 的 Elovich 方 程 拟 合 
3 em 粒 径 灰 绿 板 岩 K 元 素 的 累积 释放 量 呈 极 显 著 
相关 。3 cm 1 cm 粒 径 灰 绿 板 岩 累积 释放 量 的 4 分 
布 范围 分 别 为 : 80.065% ~91.707% 、89.520% ~ 
94.999%, 均 达到 80% 以 上 ,说 明 采 用 修正 的 Flovich 
方程 能 够 用 来 预测 两 种 粒 径 灰 绿 板 岩 矿 质 元 素 的 
实际 释放 过 程 。 相 对 于 1 cm 粒 径 灰 绿 板 岩 的 释放 
模拟 ,修正 的 了 lovich 方 程 描述 3 cm 粒 径 灰 绿 板 岩 的 
矿质 元 素 累 积 释 放 过 程 更 精准 。 

修正 的 Elovich 方 程 可 以 用 来 描述 矿质 元 素 的 
淋 溶 过 程 , 当 w=1 时 ,y=a, 可 认为 方程 中 参数 a 表示 
矿质 元 素 在 干 湿 冻 融 循环 1 次 后 的 淋 深 量 ,也 可 近 
似 认为 参数 a 为 初始 淋 溶 量 。 从 表 3 可 得 两 种 粒 径 
灰 绿 板 岩 经 过 拟 合 后 方程 中 的 参数 a 均 为 负 值 ,说 
明 经 过 1 次 干 湿 冻 融 循 环 后 灰 绿 板 岩 中 的 矿质 元 素 
发 生 了 淋 失 但 未 达到 淋 溶 量 的 累积 。 由 方程 y=a+ 
binx 可 得 矿质 元 素 淋 溶 速率 dy/dx=b/x ,说 明 是 影 
响 矿质 元 素 淋 溶 速率 的 参数 。 表 3 中 所 有 矿质 元 素 
拟 合 后 方程 中 参数 5 均 大 于 1, 当 x>1 时 ,5b 值 增 大 ， 
blac 也 增 大 ,dy/dx IL b Œ YER XS ; 24 x=1 时 ,dy/dx=0b， 
此 时 参数 5 表示 x=1 时 的 矿质 元 素 释 放 速 率 ,同时 也 
是 最 大 释放 速率 。 

FH dy/dx-b/x 可 知 ,矿质 元 素 淋 溶 速率 与 方程 中 
的 参数 5b 有关, 当 5 增 大 时 ,矿质 元 素 的 淋 溶 速率 也 
相对 较 大 ;说 明 在 同一 循环 次 数 (xz 1 时 ) ,不 同 粒 径 
矿质 元 素 的 淋 溶 速率 与 参数 有关,5 值 越 大 ,矿质 
元 素 淋 溶 速率 越 大 。 由 表 3 可 得 ,两 种 粒 径 灰 绿 板 
岩 各 矿质 元 素 淋 溶 采用 修正 的 ovich 方 程 拟 合 后 ， 


同一 元 素 不 同 粒 径 拟 合 后 的 5 值 存在 较 大 差异 , 具 
体 表 现 为 同一 元 素 粒 径 越 小 ,2 值 越 大 ,说 明 粒 径 小 
的 灰 绿 板 岩 中 矿质 元 素 的 淋 深 释放 速率 越 大 。 
2.3.3 双 常 数 速 率 方 程 (Freundlich 修正 式 ) WH 
数 速率 方程 即 是 Freundlich 修正 式 , 可 描述 能 量 分 
布 不 均匀 的 反应 过 程 ,常用 于 反应 过 程 较 复杂 的 释 
放 动 力学 过 程 ”, 通 过 双 常 数 速率 方程 模拟 不 同 粒 
径 灰 绿 板 岩 矿 质 元 素 累 积 释放 变化 ,结果 见 表 4。 

灰 绿 板 岩 各 人 矿质 元 素 的 淋 溶 释放 模拟 采用 双 
常数 速率 方程 拟 合 均 达 到 显著 性 水 平 ,其 中 两 种 粒 
径 S 元 素 的 释放 拟 合 达 到 极 显 著 水 平 ,说 明 双 常数 
速率 方程 可 更 好 的 表征 S 元 素 的 累积 释放 过 程 。3 
cm 粒 径 灰 绿 板 岩 矿 质 元 素 累积 释放 过 程 拟 合 的 忌 
分 布 范围 为 0.882~0.981 ,均值 为 0.915; 1 em 粒 径 灰 
绿 板 岩 拟 合 的 尼 分 布 范围 为 0.774~0.946 ,均值 为 
0.836, 可见 双 常数 速率 方程 对 3 cm ALARA AD 
质 元 素 的 释放 过 程 拟 合 较 好 ,可 能 是 因为 3 em fie 
灰 绿 板 岩 在 干 湿 冻 融 循环 过 程 中 与 水 接触 面积 的 
均匀 程度 较 1 cm 粒 径 灰 绿 板 岩 差 , 板 岩 与 水 作用 的 
强度 存在 差异 大 。 
2.3.4 HHZ Chute 和 Quizk 根 据 基本 扩散 作 
用 ,导出 扩散 量 与 时 间 的 关系 式 , 并 对 于 较 短 的 扩散 
过 程 经 过 简化 后 得 到 现在 常 运 用 的 抛物 线 方程 *， 
该 方程 常用 于 描述 由 多 个 扩散 机 制 控制 的 反应 过 程 ， 
最 适合 描述 物质 在 颗粒 内 部 扩散 过 程 的 动力 学 ,对 
于 颗粒 表面 物质 的 扩散 过 程 往往 不 适合 。 

不 同 粒 径 灰 绿 板 岩 矿质 元 素 释放 过 程 采用 抛 
物 线 方程 拟 合 ,结果 见 表 5。 采 用 抛物 线 方程 拟 合 
两 种 粒 径 灰 绿 板 岩 的 尼 均 大 于 0.8 , 随 循环 次 数 的 增 
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表 4 双 常 数 速率 方程 拟 合 


Tab.4 Fitting of the double constant rate equation 


粒 径 /cm 元 素 方程 Uca UU) 平均 精度 (4)/% 
3 Ca Iny=2.639+0.914Inx 0.890° 66.020 
K Iny=3.744+0.813Inx 0.926 74.912 
Mg Iny=2.964+0.828Inx 0.883° 55.968 
P 1nyz0.34640.7671nx 0.895 79.026 
S Iny=-0.225+0.9271nx 0.981” 87.140 
1 Ca lny=4.236+0.777lnx 0.815 69.131 
K Iny=5.197+0.728Inx 0.837 83.409 
Mg Iny=4.002+0.799Inx 0.774 74.505 
P Inyz1.51240.683lnx 0.808" 78.732 
S Iny=0.005+0.97 Lux 0.946" 95.116 
表 5 抛物 线 方程 拟 合 
Tab.5 Fitting of parabolic equation 
粒 径 /cm 元 素 方程 决定 系数 (RR) 平均 精度 (4)/% 
3 Ca y=-445.764+135.948x" 0.940° 75.084 
开 y=-662.617+240.866x"° 0.958" 80.564 
Mg y=-354.603+121.627x°° 0.925 69.751 
P y=-16.498+6.412x°° 0.928" 83.271 
S y=-22.843+7.642x°° 0.987" 94.771 
1 Ca y=-950.494+342.39x°° 0.861 81.689 
K y=-1672.109+683.198x°° 0.874 89.181 
Mg y=-910.177+308.029x°° 0.828 82.567 
P y2-29.697413.5834? 0.841 84.340 
S y=—40.156+12.306x” 0.968" 94.561 


加 ,3 cm PLETAL SIR VY) FRE ER 
际 值 与 预测 值 之 间 有 较 好 的 相关 关系 。1 cm HILFE 
灰 绿 板 岩 矿质 元 素 的 释放 过 程 显著 性 较 差 , 仅 S 元 
素 为 差异 显著 。1 cm 粒 径 灰 绿 板 岩 采用 抛物 线 方 
程 拟 合 4 均 高 于 80% ,两 种 粒 径 灰 绿 板 岩 中 S$ 元 素 
的 尼 和 4 很 高 ,说 明 抛 物 线 方程 可 表征 S 元 素 的 释 
放 动 力学 过 程 。 

2.3.5 矿质 元 素 释放 动力 学 方程 优选 ”不 同 矿质 元 
素 的 性 质 导 致 其 在 干 湿 冻 融 循环 过 程 中 释放 的 动 
力学 机 制 有 所 不 同 。 从 拟 合 的 4 来 评价 ,Ca Mg, K, 
P 元 素 累 积 释 放 过 程 的 方程 拟 合 效果 依次 均 为 : 修 
正 的 Elovich 方 程 > 抛物 线 方程 > 双 常 数 速率 方程 > 一 
级 动力 学 方程 。Ca KMg, .这 4 种 元 素 的 释放 模型 
与 修正 的 Elovich 方 程 拟 合 较 好 , 尼 较 高 , 且 进 行 下 
检验 发 现 函 数 拟 合 的 显著 水 平 较 高 。 通 常 Elovich 
方程 描述 活化 能 变化 较 大 的 反应 过 程 ,说 明 干 湿 冻 
融 循 环 过 程 中 这 4 种 元 素 从 灰 绿 板 岩 淋 溶 释 放 的 过 


程 中 活化 能 变化 较 大 且 相 对 较 复杂 。 双 常数 速率 
方程 和 抛物 线 方程 均 可 较 好 的 拟 合 S 元 素 的 累积 释 
放 过 程 。 

两 种 粒 径 灰 绿 板 岩 的 元 素 累 积 淋浴 释放 动力 
学 方程 拟 合 优 劣 均 呈 现 修正 的 Elovich 方程 > 抛物 线 
方程 > 双 常 数 速率 方程 > 一 级 动力 学 方程 ,修正 的 
Elovich 方 程 描 述 灰 绿 板 岩 的 矿质 元 素 释 放 更 能 
映 累 积 释 放量 随时 间 的 变化 趋势 , 即 用 修正 的 Blov- 
ich 方 程 可 以 预测 灰 绿 板 岩 在 某 一 循环 次 数 下 的 累 
只 释放 量 或 者 矿质 元 素 的 供给 量 。 


3 讨论 


许多 学 者 在 研究 岩石 风化 过 程 中 发 现 岩石 内 
部 元 素 会 发 生 一 定 程 度 的 淋浴 并 改变 土壤 元 素平 
衡 W ,但 关于 不 同 粒 径 灰 绿 板 岩 中 矿质 元 素 的 释 
放 规 律 鲜 有 报道 。 本 文通 过 研究 发 现 经 过 80 次 冻 
融 干 湿 循环 后 1 cm 粒 径 灰 绿 板 岩 矿 质 元 素 淋 失 量 
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大 于 3 cm 粒 径 ,两 种 粒 径 灰 绿 板 岩 经 历 干 湿 冻 融 循 
环 后 表面 微观 结构 如 图 4 所 示 ,1 em 粒 径 灰 绿 板 岩 
表面 的 碎 悄 颗粒 、 微 裂 际 数量 明显 多 于 3 cm ALE 
通过 PCAS 软件 对 板 岩 表面 微观 结构 分 析 计 算得 出 
1 em,3 cm 粒 径 灰 绿 板 岩 的 孔 院 率 分 别 为 28.819% 、 
22.72% ,了 筷 际 越 多 可 为 水 与 各 种 矿物 之 间 提 供 反 应 
的 机 会 越 多 , 且 本 身 比 表面 积 较 大 ,更 易于 矿物 元 
素 的 溶解 释放 ” ,使 元 素 淋 溶 量 增 大 。 

本 文通 过 对 不 同 动力 学 方程 拟 合 后 得 出 Ca、K、 
Mg、P 元 素 的 累积 释放 最 优 动 力学 方程 为 修正 的 
Elovich 方 程 。 其 中 修正 的 Elovich 方 程 最 能 体现 久 元 
素 的 累积 释放 过 程 , 尼 达到 0.93 以 上 ,与 王 瑾 等 
的 研究 结果 大 致 相同 。 修 正 的 Elovich 方程 能 较 好 
地 描述 元 素 的 累积 释放 过 程 可 能 是 由 于 灰 绿 板 岩 
的 矿物 成 分 中 云母 和 长 石 含量 相对 较 多 ,云母 类 矿 
物 一 般 为 2:1 型 层 状 结构 ,对 层 间 K 元 素 的 束缚 能 力 
较 弱 ,使 K 元 素 与 离子 半径 和 水 化 能 较 大 的 Ca”* 进 行 
交换 ,或 被 Na 置换 ”等 增加 反应 过 程 的 复杂 性 ” ; 
Ca 相对 活泼 , 干 湿 冻 融 循环 过 程 中 Ca 元 素 既 可 能 
发 生 在 斜 长 石 中 与 Na 元 素 的 类 质 同 像 置换 ,还 会 以 
碳酸 盐 的 形式 溶解 于 水 中 等 复杂 反应 ;在 灰 绿 板 
岩 中 Mg 元 素 主要 来 源 于 绿 泥 石 中 ,通常 以 与 CI 化 
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合 反应 或 与 SOT 相 结 合 的 方式 溶解 释放 ,在 碱 性 环 
境 中 镁 碳酸 盐 与 钙 碳酸 盐 还 会 沉 诈 形成 钙 与 镁 的 
复 盐 ”; 冻 融 干 湿 循环 过 程 中 P 元 素 的 行为 与 岩石 
中 磷 矿 物 的 种 类 和 这 些 矿物 在 风化 溶解 时 的 溶解 
度 有 关 , 同 时 灰 绿 板 岩 中 粘土 矿物 对 P 的 吸附 作用 
还 影响 着 冻 融 循环 过 程 中 P 元 素 的 释放 浓度 "。 由 
此 可 见 ,Ca、K、Mg、P 元 素 在 冻 融 干 湿 循 环 过 程 中 的 
淋 溶 释放 是 个 复杂 的 反应 过 程 ,因此 可 采用 修正 的 
Elovich 方 程 来 描述 各 矿质 元 素 的 释放 动力 学 过 程 。 

采用 修正 的 Elovich 方 程 可 以 预测 不 同 粒 径 板 
岩 在 某 一 循环 次 数 下 的 矿质 元 素 释 放量 ,有 助 于 预 
测 不 同 循环 次 数 ( 压 砂 年 限 ) 的 矿质 元 素 供给 量 ,再 
结合 不 同 作 物 不 同 生育 期 对 各 矿质 元 素 的 需求 量 ， 
便 可 通过 调节 灌水 频次 等 来 满足 作物 所 需 矿 质 元 
素 的 含量 ,做 到 改善 压 砂 地 土壤 矿质 元 素平 衡 的 同 
时 为 作物 提供 合适 的 矿质 元 素 ,做 到 压 砂 地 土壤 矿 
质 元 素 丰 缺 调控 与 作物 营养 调控 相 结合 。 


4 结 论 


-H 


(1) 干 湿 冻 融 作用 下 灰 绿 板 岩 矿质 元 素 的 释放 
有 快速 反应 、 趋 于 稳定 释放 2 个 阶段 ,其 中 快速 反应 


图 4 两 种 粒 径 灰 绿 板 岩 表面 微观 形 貌 图 (1000 售 ) 


Fig.4 Surface microscopic morphology of gray-green slate with two particle sizes (1000 times) 
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阶段 1 em 粒 径 灰 绿 板 岩 的 释放 速 
灰 绿 板 岩 的 释放 速率 。 

(2) 灰 绿 板 岩 粒 径 越 小 ,颗粒 比 表 面积 越 大 ,与 
水 接触 越 充 分 ,岩石 表明 可 溶性 元 素 易 流失 ,导致 
其 矿质 元 素 的 淋 溶 释放 量 越 大 。 因 此 ,在 调控 压 砂 
地 土壤 元 素平 衡 时 可 适当 考虑 覆 砂 砾石 的 粒 径 大 
小 分 布 。 

(3) 修正 的 Elovich 方 程 能 更 好 地 描述 灰 绿 板 
岩 的 累积 释放 特征 ,其 尺 和 4 分 别 为 :0.895~0.988 
和 80.065%~94.999%, 均 显著 相关 ,说 明 干 湿 冻 融 作 
用 下 矿质 元 素 的 释放 过 程 是 由 多 个 因素 控制 的 复 
杂 反 应 过 程 。 


束 率 高 于 3 cm 粒 径 
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Mineral element release characteristics and release models for gray-green slate 
under alternating freeze-thaw and dry-wet conditions 


WANG Jie, LI Wangcheng"", MU Min', DONG Yaping 
(1. School of Civil and Hydraulic Engineering, Ningxia University, Yinchuan 750021, Ningxia, China; 2. Engineering Research 
Center for Efficient Utilization of Modern Agricultural Water Resources in Arid Regions, Ministry of Education, Yinchuan 750021, 
Ningxia, China; 3. State Key Laboratory of Land Degradation and Ecological Restoration in Northwest China, Yinchuan 750021, 


Ningxia, China) 


Abstract: In response to the problems of soil fertility decline in gravel-mulched land in the arid zone of central 
Ningxia Hui Autonomons Region, China and unclear leaching and release rules of gray- green slate mineral 
elements, this study conducted an indoor simulated freeze-thaw and dry-wet cycle experiment of gray-green slate, 
used the modified Elovich equation, parabolic equation, double constant rate equation, and first- order kinetic 
equation to fit the cumulative release curve of each mineral element, and study the release kinetic characteristics 
and optimal fitting kinetic equation of slate with different particle sizes under the action of dry-wet freeze-thaw 
cycles. The results show that the total leaching of mineral elements in the leaching solution of the gray- green 
slates of two-grain sizes increases with the number of cycles, and the cumulative total leaching of mineral 
elements and leaching rate of 1 cm grain size slate are greater than those of 3 cm grain size slate at different cycle 
times. The release of mineral elements from the gray-green slate was a physical and chemical process controlled 
by multiple factors, and the release process could be roughly divided into the rapid reaction stage and the reaction 
to the equilibrium stage. The supply of Ca, K, Mg, and P elements in the gravel-mulched land soil was best predicted 
by the modified Elovich equation, and the parabolic equation was most suitable for describing the release pattern 
of S elements. The findings of this study could serve as a decision reference for soil fertility regulation of dry 
farmlands in the arid zone of central Ningxia. 
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